Hot phonon decay in supported and suspended exfoliated graphene 
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Near infrared pump-probe spectroscopy has been used to measure the ultrafast dynamics of pho- 
toexcited charge carriers in monolayer and multilayer graphene. We observe two decay processes 
occurring on 100 fs and 2 ps timescales. The first is attributed to the rapid electron-phonon ther- 
malisation in the system. The second timescale is found to be due to the slow decay of hot phonons. 
Using a simple theoretical model we calculate the hot phonon decay rate and show that it is signifi- 
cantly faster in monolayer flakes than in multilayer ones. In contrast to recent claims, we show that 
this enhanced decay rate is not due to the coupling to substrate phonons, since we have also seen 
the same effect in suspended flakes. Possible intrinsic decay mechanisms that could cause such an 
effect are discussed. 
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The symmetric, linear electronic band structure of 
graphene gives rise to some very unusual physical prop- 
erties, such as quantised transmission^, extremely high 
thermal conductivity^ and high carrier mobility^ An im- 
portant underlying feature is the very strong electron- 
phonon coupling that exists in graphene, which is re- 
vealed by the presence of Kohn anomalies^ In graphite, it 
is known that strongly coupled optical phonons have high 
quantum energies of up to 0.2 eV and are excited only 
by electrons of elevated energy4^ To progress towards 
applications in real (high-current) circuits and devices, 
it is crucial to understand how graphene behaves under 
such high energy, non-equilibrium conditions. Despite 
the surge of interest in this material and its potential 
applications, investigations into the kinetic properties of 
"hot" charge carriers remain rather limited 

Hot electron relaxation in large area, epitaxially 
grown graphene layers using pump-probe spectroscopy 
has been studied previously^— These measurements 
point to biexponcntial decay dynamics characterised by 
a fast ~100fs component and a slow ~2ps component. 
There is, however, significant variation in the reported 
timescales. Epitaxial graphene exhibits inhomogeneity in 
layer thickness on the micron scaled and can result in sig- 
nificant variability in relaxation dynamics from sample to 
sample. 8 Pump-probe measurements have also been per- 
formed on mechanically exfoliated graphene ) 17 i 18 which 
is homogeneous over much greater length scales. It was 
concluded that the slow relaxation process was caused by 
the coupling to phonons in the substrate*^ 

In this paper we use near infrared pump-probe spec- 
troscopy to investigate the relaxation dynamics of hot 
carriers in mechanically exfoliated graphene. Similar to 
previous results we find that the relaxation occurs on two 
timescales, one fast (~100fs) and the other slow (~2ps). 
By measuring the relaxation in monolayer and multi- 
layer graphene flakes we show a clear correlation between 
the slow decay rate and flake thickness, with the fastest 
rate observed for monolayer graphene. This slow decay 
rate is found to occur in both supported and suspended 
flakes. Therefore, in contrast to recent claimS ) 17 i 19 the 
enhanced relaxation is not due to the coupling to sub- 
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FIG. 1. (Color online) Differential reflectivity (nomalised at 
t — 0) for a 10-layer sample as a function of the probe delay at 
different probe energies. Inset: the electronic band structure 
of graphene at the K point of the Brillouin zone with the 
relevant pump and probe photon energies shown. 

Samples were prepared from mechanically exfoliated 
natural graphite deposited on 100 /im thick glass sub- 
strates. The number of layers in each sample was de- 
termined by optical contrast^ and four-wave mixing 
measurements <2i Raman spectroscopy was used to con- 
firm the thickness of monolayers and bilayers from the 
profile of the 2D peak and the ratio of the G and 2D 
peak intensities^ Suspended samples were made by de- 
positing flakes on glass substrates that had an array of 
holes etched into the surface. Square holes between 1.5 
and 4.5 /zm were fabricated by reactive ion etching to a 
depth of 290 nm. 

To perform the pump-probe measurements, a 
Thsapphirc mode-locked laser, center wavelength 830 nm, 
pulse width ^180 fs, was used to optically pump the 
samples. A probe pulse, ^240 fs, was generated using 
an Optical Parametric Oscillator with variable central 
wavelengths spanning 1100- 1400 nm. The pump and 
probe pulses are set with an intensity ratio of > 10:1. 
The beams were aligned into a microscope and raster 
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scanned through a 1.2 NA water immersion lens onto the 
sample, ensuring a spot size < 1.5 /mi, which is signif- 
icantly smaller than the graphenc flakes. We estimate 
the intensity of the combined beams at the sample to be 
~ 0.03 J/m . Transient dynamics of the photoelectrons 
are obtained by measuring the reflected probe light, fo- 
cussed through 1100 nm long-pass filters onto a large area 
photodiode, and recorded as a function of delay between 
pump and probe signals. The change in the probe light 
reflectivity is observed to be linearly dependent upon the 
pump intensity. By varying the focus of the microscope 
objective lens, we find that the dynamics measured are 
independent of the size of the excitation region, indicat- 
ing that lateral charge diffusion does not affect our re- 
sults. 

Figure [T] shows a typical pump-probe measurement of 
(normalised) differential reflectivity, SR, for a 10-laycr 
thick sample. When the pump and probe pulses are 
overlapped in time the pump pulse excites electrons and 
blocks available transitions for the probe, resulting in a 
decrease in the reflection from the sample. The relaxation 
of excited carriers is then observed as a function of time 
by delaying the probe pulse with respect to the pump 
and measuring SR. These data can be described by two 
exponentially decaying curves, one with a characteristic 
rate much faster than the other. The fast timescale is 
not fully resolved in our measurement as it is limited by 
the probe width. Our interest, however, lies exclusively 
in the slow timescale which is easily resolved. 
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FIG. 2. (Color online) (a) Differential reflectivity as a func- 
tion of the probe delay for samples of different thickness, mea- 
sured using a probe wavelength of 1120 nm. (b) Extracted 
phonon lifetimes as a function of layer number. Averaged 
values per layer are also plotted with the standard deviation 
as an error bar. The dotted line is a guide to the eye. 

Wc have performed measurements as a function of the 



number of graphenc layers in a large number of different 
flakes. Figure EJa) contrasts the reflectivity measured 
in flakes from single to six layers supported on a glass 
substrate. We observe a clear correlation between layer 
number and decay rate: faster relaxation occurs for the 
thinner samples. However, for ~4 layers and above we 
find that all data collapse onto a single curve. This indi- 
cates that at these thicknesses the sample behaves as a 
bulk material and the rate of relaxation is constant. This 
suggests that the mechanism that removes energy from 
the system in thin layers is not available to the thicker 
layers. 

The fast and slow relaxation timescales have been 
previously attributed to electron-electron and electron- 
phonon interactions, respectively^ However, it is well 
known that clcctron-phonon coupling in graphene 
features Kohn anomalies^ that lead to very fast 
thermalisation with the phonon bath.— Indeed, re- 
cent calculation ;) 7 ' 23 indicate that thermalisation with 
phonons is completed within the first ~T00 fs following 
photoexcitation. This suggests that the longer timescale, 
of the order picoseconds, is associated with the relaxation 
of the thermalised electron-phonon bath through the de- 
cay of hot phonons. 

In order to better understand the mechanisms of the 
decay, we compare our data to a simple coupled rate 
equation model which estimates the electron and optical 
phonon temperatures after photoexcitation from ultra- 
fast optical and near-infrared pulses. This model is based 
on earlier ones developed in Refs. [lj| and [HI which we 
outline briefly below. We calculate the phonon genera- 
tion rate using the relation 

/oo 
E(E-hu> ph ) x [ Pe -p a ]dE, (1) 
-oo 

where w p h is the optical phonon frequency (^180meV), 
p c = f(E)(l — f(E — fiujp} 1 ))(npi l + l) is the probability of 
emitting a phonon, p a = f(E~huip^)(l — f(E))n p ^ is the 
probability of absorbing a phonon and /(E) is the Fermi- 
Dirac distribution of electrons for a given electron tem- 
perature, T i. Here a = 9/2 x j3 2 \irpujp\fi , is the 
electron-phonon coupling strength, where p is the density 
of graphene (~ 7.6 x 10~ 7 kg/m 2 ), vp is the Fermi veloc- 
ity and j3 = 45eV/nm^ The phonon occupation num- 
ber, n p h, is related to the electron temperature through 
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where C e \ is the electronic heat capacity— (Note that 
the factor 2 arises from the presence of both T and K 
point phonons. Here we treat these phonons as being 
equivalent due to their similar frequencies.) The rate of 
optical phonon generation is governed by the equation 
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where n ph is the phonon occupation at room temper- 
ature, To, and T p h is the phonon decay rate. The 
second term in Eq. |3] describes the anharmonic de- 
cay of hot phonons. M p h is the number of T and K 
point phonons (per unit area) that are able to couple 
to the hot electrons, which is estimated by consider- 
ing the maximum and minimum momentum possible for 
an emitted T or K point phononJ^ This gives M p h = 
(2/4?r) [(V2E ma ^/hv F ) 2 - (wph/wF) 2 ] , where E max is up- 
per energy of the hot electons that are able to emit 
phonons, found by calculating T p h for a given T c \. (Note 
that the factor of V2 in front of E max is due to the con- 
servation of pscudospin which limits scattering to states 
within the same semi-cone in the Brillouin zone, and the 
factor 2 arises from valley degeneracy.) 

Using Eq. [2] we can evaluate the temperature of the 
electron bath as a function of time after excitation. Since 
the energy from our excitation pulse is dissipated through 
intraband optical transitions, an initial electron temper- 
ature T o i(0) is calculated from the absorbed fluence and 
Coi- We account for cooling during the excitation pulse 
by estimating an average emission of 5 optical phonons 
per excited electron during the excitation pulsed Given 
these starting conditions, a time dependent reflection 
change is found using 

M(t)KtMl (_i_)_ tanh (_|y , (4) 

which gives the change in the strength of interband tran- 
sitions of the probe pulse due to the electron temperature 
and the probe photon energy, E pr . We use probe ener- 
gies ranging from 0.94- 1.11 eV (1120 - 1320nm). After 
convolving Eq. @]with the correlation of our pulses, and 
using just the absorbed pump flucncc and r p h as fitting 
parameters, we are able to fully describe the time depen- 
dent change in reflection measured in our experiments. 

The two data sets shown in Fig. [T] are measured with 
two different probe energies on the same flake: 0.94 eV 
and 1.11 eV. The fast initial decay is roughly the same 
for both energies (< 200 fs), whereas the slower decay de- 
pends on the probe energy. The fits (solid lines) give 
T p h=3.35 ±0.1 ps and T p h=3.6 ±0.05 ps for probe energies 
of 0.94 eV and 1.11 eV, respectively. This slow relaxation 
corresponds to loss of energy from the hot phonon pop- 
ulation, presumably through anharmonic decay. This 
creates a bottleneck for energy loss. Without this fi- 
nite decay mechanism, hot phonons would accumulate to 
high occupation numbers, thus changing the timescalcs 
of electronic kinetic processes and leading to low damage 
thresholds in, for example, high current graphene devices. 

One can clearly see different dynamics measured for 
different probe energies. This behaviour arises within 
our model due to the different electron energies probed in 
each measurement: higher probe energies, probing higher 
energy electrons, will yield faster decay dynamics. Our 
model, meanwhile, yields similar values for the phonon 
decay time. Theoretical calculations of the hot optical 
phonon lifetimes in graphene due to anharmonic decay 



into acoustic phonons were carried out in Ref . [26l and val- 
ues between 2 and 4 ps were reported for both zone-center 
E 2g and zone-edge A\ modes for phonon temperatures in 
the 500 — 900 K range. These calculations are clearly con- 
sistent with our measurements. Moreover, the values we 
obtain are similar to those obtained from recent pump- 
probe measurements on epitaxially grown graphene using 
a similar models 

The change in probe reflection as a function of pump 
delay for our supported flakes of different thickness is 
shown in Fig. HJa). Similarly to Ref. [l9j, the extracted 
values of T p h from fits to the above theory (solid lines) in- 
crease as a function of increasing layer number, Fig.^b). 
The authors suggested that this is a result of the coupling 
between charge carriers in graphene and phonons in the 
underlying substrate. In order to investigate the role of 
the substrate in this effect, the next logical step is to 
compare the relaxation behaviour observed in supported 
to that in suspended graphene flakes. 

Figure EIa,b) show results obtained from measure- 
ments on flakes suspended over 3.5 (im holes. Flakes 
were characterised by atomic force microscopy to estab- 
lish that they were truly suspended over the holes in the 
substrates, Fig. (3{c). The interaction with the substrate 
phonons will decay on a lengthscale comparable to the 
phonon wavelength. The maximum wavelength possi- 
ble for an emitted substrate phonon due to coupling to 
electrons in graphene is (2-nvp / w' ph ) , where w ph is the 
substrate optical phonon frequency. For normal opti- 
cal phonon energies (^100 mcV) the interaction typically 
decays over distances < 100 nm. Since the depth of the 
etched holes in our case is 290 nm, we expect interactions 
with substrate phonons to be substantially suppressed as 
a relaxation route. In our pump-probe measurements, 
the high spatial resolution (<1.5/im) also guarantees 
that only the suspended portions of the flakes are in- 
vestigated. 

Time resolved measurements of the suspended flakes 
are shown in Fig. EJa). Flakes of different thicknesses 
were measured under the same experimental conditions 
as the supported flakes. We clearly observe a similar layer 
dependence as observed for supported layers, with faster 
decay times observed for thinner flakes. The extracted 
phonon lifetimes are also similar to those of the supported 
flakes. In Fig. [2Kb) we compare the suspended measure- 
ments to the average value of the supported flakes (with 
the standard deviation as error bars). Indeed, it would 
appear that relaxation of suspended monolayer graphene 
may be even faster than that for the supported flakes. 
This strongly suggests that the substrate is not respon- 
sible for the fast decay observed in thin flakes, and that 
the mechanism of phonon relaxation in graphene has an 
intrinsic origin. 

A possible explanation lies with out-of-plane (flexural) 
phonons, which are suppressed in thicker layers. One 
would expect flexural modes to be present for monolayer 
graphene supported on a rough glass substrate as the 
phonon wavelengths in question are comparable to the 
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FIG. 3. (Color online) (a) Differential reflectivity for sus- 
pended mono- and few-layer graphene. (b) Extracted phonon 
lifetimes for mono- and few-layer suspended graphene plotted 
alongside the averaged values extracted from the supported 
flakes, (c) An AFM image of a suspended flake covering a 
2.5 pm, 290 nm deep hole. 

roughness of the substrate. Out-of-planc modes are sup- 
pressed in thicker layers due to interactions between the 
graphene planes. As the number of layers increases one 
would expect an asymptote in the observed dynamics as 
the system begins to resemble a rigid body. From our 



data, we suggest that this asymptote is already reached 
at approximately four layers. 

In conclusion, near infrared pump-probe spectroscopy 
has been used to measure the ultrafast dynamics of pho- 
toexcitcd charge carriers in monolayer and multilayer 
graphene. We observe two timescales in the decay oc- 
curring on ~100fs and ~ 2 — 3ps. The fast relaxation 
timecale can be attributed to electron-phonon thermal- 
ization, whereas the slower timescale represents a bottle- 
neck in the relaxation process due to anharmonic decay 
of hot phonons. Using a simple theoretical model we have 
calculated the decay rates in different graphene flakes and 
shown that the hot phonon decay is faster in monolayer 
than in multilayer graphene. Comparing our results on 
supported and suspended flakes we have demonstrated 
that substrate phonons are not the mechanism for re- 
moving energy from the system. The possible origin of 
the intrinsic mechanism is the enhanced coupling to out- 
of-planc flcxural phonon modes. 
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